%HCR signal isolation comparison between old method using Fiji and new
%Matlab method.

%Step1 load arrays of same slice processed with old and new method. 1
%channel at a time.

figure;histogram(old(new>100),'Normalization','cdf','DisplayStyle','stairs')
hold on
hold on;histogram(new(new>100),'Normalization','cdf','DisplayStyle','stairs')
xlim([0 5e3])
xticks([0 0.5 1])
xticks auto
xlim([0 5e3])
yticks([0 0.5 1])
set(gca, 'TickDir', 'out')
xticks([0:1e3:5e3])
box off